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PREFACE 

This  report  describes  physical  and  mathematical  models  which  form  the 
basis  for  two  digital  computer  codes  which  can  be  used  to  simulate  the 
motion  of  a Short  Range  Air  Defense  System  (SHORADS)  rocket  before  it  is 
acquired  by  a launcher  fixed  tracking  radar.  Listings  of  the  two  computer 
codes  as  well  as  typical  simulation  results  are  also  included  herein. 

Much  of  the  material  in  this  report  also  constitutes  part  of  Mr.  James 
R.  Beaty's  Master  of  Science  thesis. 

The  assistance  of  Mr.  John  Howerton  of  the  Army  Missile  Research  and 
Development  Command,  Redstone  Arsenal,  Alabama,  who  provided  technical  data 
for  this  study  is  gratefully  acknowledged. 

Dr,  Joe  W.  Reece,  Co-Project  Leader,  is  presently  on  leave  from 
Auburn  University  and  did  not  participate  in  the  writing  of  this  report. 
However,  his  contributions  to  the  overall  effort  during  its  initial  stages 
were  indispensable  in  regard  to  the  successful  completion  of  this  study. 

The  authors  also  gratefully  acknowledge  the  patience  and  expertise  of 
Mrs.  Marjorie  McGee,  who  typed  the  manuscript. 
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SECTION  1.  INTRODUCTION 
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1 . 1 General  Comments 

The  rising  costs  of  design,  development  and  testing  of  modern  weapon 
systems  coupled  with  the  availability  of  high-speed  digital  computers  have 
led  to  increased  utilization  of  digital  computer  codes  to  simulate  proposed 
and  existing  weapon  systems.  A wealth  of  information  about  the  operation 
of  a weapon  system  can  be  obtained  from  such  simulation  codes  if  they  are 
properly  devised.  The  type  of  information  may  range  from  results  which  are 
only  order-of-magnitude  estimates  of  say,  accelerations  and  angular  rates, 

\ 

to  very  accurate  predictions  or  reproductions  of  the  system  behavior  as  a 
function  of  time.  This  report  describes  the  development  and  use  of  two 
simulation  codes  of  moderate  complexity  which  are  intended  for  use  in  de- 
termining the  effects  of  various  anomalies,  such  as  dynamic  imbalance  and 

J 

winds,  and  of  certain  inherent  characteristics  of  a particular  type  of 
weapon  system. 

1.2  Scope  of  the  Present  Study 

This  report  deals  with,  the  problem  of  simulating  the  motion  of  a Short 
Range  Air  Defense  System  (SHORADS)  type  rocket  during  the  "capture"  or  pre- 
radar-acquisition portion  of  flight.  This  portion  of  the  rocket's  flight 
includes  its  launch  from  a tube  attached  to  the  rotating  turret  of  a ground- 
based  mobile  launcher  and  the  rocket's  motion  prior  to  its  being  acquired  by 
a launcher-based  tracking  radar.  The  part  of  the  flight  after  launch  and 
prior  to  radar  acquisition  can  be  divided  into  two  subparts.  The  first 
subpart  is  a stabilized,  but  unguided,  flight  phase  which  begins  when  the 
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aft  end  of  the  rocket  leaves  the  launch  tube  and  ends  when  the  aft  end  of 
the  rocket  enters  the  field  of  view  (FOV)  of  an  infrared  (IR)  angle  sensor, 
the  sight  axis  of  which  is  bore-sighted  with  that  of  the  radar.  (The 
IR  sensor  senses  the  position,  within  its  conical  FQV,  of  flares  on  the 
rocket's  aft  end.)  The  second  subpart  is  a flight  phase  in  which  signals 
from  the  IR  sensor  and  pre-specif ied  range  versus  time  relationships  are 
used  to  determine  guidance  commands  for  the  rocket.  This  phase  begins  when 
the  IR  sensor  acquires  the  rocket's  flares  and  ends  when  the  tracking  radar 
acquires  contact  with  the  rocket. 

A simulation  code  has  previously  been  developed1  for  the  purpose  of 
simulating  the  motion  of  a SHORADS  type  rocket  from  launch  tube  exit  until 
the  end  of  its  flight.  That  code  is  very  complicated,  due  to  the  fact  that 
an  attempt  was  made  to  simulate  the  actual  operation  of  all  subsystems.  It 
does  not,  however,  contain  a provision  for  simulating  the  launch  phase,  nor 
does  it  contain  a rocket  mathematical  model  which  includes  certain  anomalies, 
considered  in  this  study  and  included  in  the  "post-launch"  simulation  code 
described  herein. 

Two  digital  computer  codes  were  developed  during  this  study.  The  first 
is  a simplified  launch  dynamics  code  and  the  second  i s the  post-launch 
simulation  code  cited  above. 

Section  2 of  this  report  contains  a description  of  the  physical  and 
mathematical  models  which  form  the  basis  for  the  launch  dynamics  code.  The 
launcher  is  modeled  as  a rigid,  fixed  base  to  which  a rigid  turret  with  two 
rigid  launch  tubes  is  attached.  The  turret  is  assumed  to  rotate  about  a 
vertical  axis  with  constant  angular  speed.  For  launch  dynamics  purposes, 
the  rocket  is  modeled  as  a constant  mass,  axisymmetric  rigid  body.  The 
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effects  of  tipoff  of  the  rocket  are  modeled.  The  desired  output  of  the 
launch  dynamics  code  is  the  initial  state  vector  of  the  rocket  at  the 
instant  its  aft  end  leaves  the  launch  tube. 

The  physical  and  mathematical  models  of  the  rocket  are  presented  in 
Section  3.  The  rocket  is  modeled  as  a variable  mass  body  in  the  sense 
that  the  indirect  effects  of  variations  in  its  mass  and  moments  of  inertia 
with  time  and  the  direct  effects  of  thrust,  torque  due  to  thrust,  and  jet 
damping  due  to  internal  flow  are  modeled.  The  rocket  modeled  herein  is  of 
the  folding-wing,  cruciform,  boost-sustaln  type.  It  is  assumed  to  be  spun- 
up  to  a relatively  low  roll  rate  (say,  5 cps  ) during  launch  and  to  maintain 
an  essentially  constant  roll  rate  after  launch.  The  thrust  of  the  rocket 
as  a function  of  time  is  assumed  to  be  available  in  tabular  form.  Its  mass 
is  modeled  as  a piecewise  linear  function  of  time.  Also,  the  aerodynamics 
characteristics  of  the  rocket  are  assumed  to  be  available  in  tabular  form. 

For  example,  C^,  the  normal  force  coefficient,  is  assumed  to  be  available 
as  a tabulated  function  of  the  flight  Mach  number  and  the  incidence  angle  of 
the  rocket . 

Control  is  accomplished  by  deflecting  the  thrust  vector  of  the  sustainer 
rocket  motor  with  small  "control  vanes"  which  are  located  at  the  exit  plane 
of  the  sustainer  nozzle.  Control  torques  about  a single  transverse  axis 
fixed  in  the  rocket  are  generated  by  the  deflected  thrust.  Since  the  rocket 
is  rolling  rapidly,  control  torques  can  be  applied  alternately  to  control 
pitch  and  yaw  plane  motions  of  the  rocket.  The  magnitude  of  the  control 
torque  produced  by  deflecting  the  control  vanes  is  assumed  constant  (a  delay 
time  is,  however,  modeled)  and  the  "amount"  of  control  is  determined  by  the 
length  of  time  the  vanes  are  deflected  in  a given  direction;  i.e.,  the  dwell 
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A simple  control  law  is  described  in  Section  4.  Under  this  law,  the 
dwell  time  is  determined  as  a linear  combination  of  deviations  of  the 
rocket's  center  of  mass  from  the  launcher-to-target  line-of-sight , the  time 
rates  of  change  of  these  deviations  and  the  pitch  rate  of  the  rocket.  The 
choice  of  control  gains  which  appear  in  the  control  law  is  also  discussed 
in  Section  4. 

Results  obtained  using  the  launch  dynamics  (LD)  and  pre-radar-acquisition 
(PRA)  codes  are  presented  in  Section  5.  Plots  of  the  flight  paths  of  the 
rocket  for  when  its  motion  is  affected  by  various  anomalistic  factors  are 
shown  and  discussed. 

Section  6 contains  a summary  of  the  main  findings  of  this  study  and 
also  contains  some  recommendations  for  further  work  in  simulating  the  launch 
of  SHOFADS  rockets. 

Throughout  this  report,  the  words  rocket  and  missile  will  be  used  inter- 
changeably . 


SECTION  2.  LAUNCH  DYNAMICS 


2.1  Introductory  Comments 

Often,  the  launch  phase  is  neglected  when  simulations  of  rocket 
systems  are  conducted.  The  reason  for  this  neglect  often  is  a belief  on 
the  part  of  the  one  conducting  the  simulations  that  the  effects  of  the 
rocket’s  motion  during  the  launch  phase  may  be  accounted  for  by  varying 
the  initial  conditions  for  the  post-launch  motion.  This  may  be  the  case, 
but  one  should  not  be  led  to  believe  that  it  is  always  so.  The  main  fault 
of  the  process  of  simply  varying  initial  conditions  to  simulate  "bad  launches," 
tipoff,  etc.,  is  that  the  initial  conditions  may  not  be  consistent  and  hence 
may  lead  to  spurious  results. 

The  launch  phase  is  not  as  adequately  treated  herein  as  it  could  have 
been  if  the  majority  of  the  time  available  had  not  been  spent  on  the  post- 
launch phase.  At  the  beginning  of  this  study,  it  was  anticipated  that  an 
existing  launch  dynamics  code2  would  be  used  to  simulate  the  launch  phase. 

That  code  did  not  appear  suitable,  however,  since  it  does  not  contain  a 
model  for  a launcher  with  a rapidly  rotating  turret.  Rather  than  modify  it, 
a relatively  simple  physical  model  for  the  launcher/rocket  system  during 
launch  was  formulated,  the  associated  mathematical  model  derived  and  a new 
code  incorporating  this  new  model  written. 

2.2  Launcher/Rocket  System  Physical  Model 

The  physical  model  for  the  launcher/rocket  system  during  the  launch  phase 
is  shown  schematically  in  Fig.  2-1.  It  consists  of  a rigid,  fixed  base,  a 
rigid  turret  and  two  rigid  launch  tubes.  The  turret  rotates  relative  to  the 
base  about  a vertical  axis.  The  launch  tubes  may  be  rotated  to  different 
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elevation  settings.  It  is  assumed  that  during  launch  the  turret  rotates  at 
a constant  rate  and  that  the  launch  tubes  are  fixed  with  respect  to  the 
turret.  The  rocket  model  is  a rigid  (constant  mass)  axisymmetric.  body. 


Figure  2-1.  Launcher  Model 

2.3  Launcher/Rocket  System  Mathematical  Model 

The  launch  phase  is  divided  into  two  parts  for  the  purposes  of  deriving 
equations  of  motion  for  the  rocket  (No  equations  of  motion  are  needed  for 
the  launcher,  since  its  motion  is  specified  a priori.).  The  first  part 
of  the  launch  phase  is  one  in  which  the  rocket  translates  along  a tube, while 
it  spins  within  the  tube.  The  second  part  is  one  in  which  the  rocket  is 
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supported  only  at  a point  near  its  aft  end  and  hence  "tips  off";  i.e., 
rotates  relative  to  the  tube  about  a transverse  axis. 

Motion  During  the  First  Part  of  the  Launch  Phase 

During  the  constrained  translational  part  of  the  launch  phase,  the 


position  of  the  rocket's  center  of  mass  C is  given  by 

R_  = R.  + r , 

-C  -A  -C/A 


(2-1) 


where  is  a vector  from  point  T,  the  origin  of  the  turret-fixed  Tx^y^z^, 
system,  to  point  A,  the  origin  of  the  tube-fixed  Ax^y^z^.  coordinate  system 
and  rG/^  is  a vector  from  A to  C,  the  center  of  mass  of  the  rocket.  If  it 
is  assumed  that  the  turret  rotates  with  a constant  angular  velocity  w.jr'w.pk^,, 
where  k^is  a unit  vector  fixed  to  the  z^-axis,  the  acceleration  of  C is 

5c  = “tx[^tx(V  £c/a)] 


+ 2 Sc/A  + -C/A’ 


(2-2) 


where  the  small  circle  denotes  differentiation  of  the  components  of  a vector 
expressed  in  terms  of  the  tube-fixed  unit  vectors  It,  j and  kt . Since  the 
zf-  and  yt-motions  of  the  rocket  are  constrained  to  be  zero  during  this  part 
of  the  launch  phase,  only  the  x -component  of  R„  is  needed.  By  determining 
this  component  and  also  the  component  of  the  gravitational  force,  we  find 
that  the  equation  for  the  rocket's  translation  within  the  tube  is 


xt  = - g sin  <?G  + xt  uj2  cos2  $G  + FT/m 


(2-3) 


where  g is  the  acceleration  of  gravity,  $G  is  the  elevation  angle  of  the 
tube,  F ^ is  the  magnitude  of  the  thrust  force  (friction  is  neglected)  and 


m is  the  mass  of  the  rocket. 
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k 
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t 


r 


I 
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The  rocket  also  spins  within  the  tube.  We  let  p denote  the  spin, or 
roll,rate  of  the  rocket  relative  to  the  tube,  Tx  denote  the  x^-component 

(also  the  x-component,  where  x is  the  longitudinal  axis  of  the  rocket)  of 

the  torque  on  the  rocket  and  I denotes  the  rocket's  moment  of  inertia  about 

x 

the  xt~axis.  The  equation  for  rotational  motion  of  the  rocket  relative  to 
the  tube  is  then  simply 

i - Tx/Ix  . (2-4) 

The  roll  angle  of  the  rocket  relative  to  the  tube  may  be  found  from  the 
kinematic  equation 

A*  = P . (2-5) 

Tipoff  Equations  of  Motion 

For  the  purposes  of  this  study  it  is  assumed  that, when  the  rocket  has 
traveled  a specified  distance  down  the  tube,  the  transverse  rotation  of  the 
rocket  is  no  longer  constrained  by  the  tube  and  any  internal  guidance 
mechanism,  such  as  a carriage.  During  the  interval  of  time  beginning  when 
transverse  rotation  of  the  rocket  is  not  constrained  and  ending  when  the 
aft  end  of  the  rocket  exits  the  tube,  "tipoff"  occurs.  In  formulating  the 
equations  of  motion  for  this  part  of  the  launch  phase,  it  is  assumed  that 
the  rocket  may  rotate  transversely  about  a point  a on  the  circumference  of 
the  rocket,  which  point  a is  located  at  the  end  of  the  tube  (see  Fig.  2-2), 
We  let  Aip  and  AS  denote  Euler  angles  which  along  with  A<P  define  the 
relative  orientation  of  the  rocket-fixed  Cxyz  system  with  respect  to  the 
tube-fixed  Axtytzt  system,  p,  q and  r denote  the  rocket-fixed  component* of 
the  relative  angular  velocity,  the  rocket  with  respect  to  the  tube- 

. . • ■■  --  -■  mm  a — * ~ 


Figure  2-2.  Geoae try  During  Tipoff. 

fixed  system  and  1^  denote  the  rocket's  inertia  matrix.  Furthermore,  we 

define  d as  the  diameter  of  the  rocket  at  the  point  a and  assume  that  d is 

constant.  We  also  let  denote  the  distance  from  point  A to  the  fore  end  of 

E 

the  tube.  Finally,  we  define  the  following  transformation  matrices  (assuming 
A9  and  Ai|i,  but  not  A<{>,  are  small  angles): 


1 

Ai/> 

-A  6 

A = 

=r/t 

-A^  cos  A<f>  + A0  sin  A4> 

cos  A(j) 

sin 

A<(> 

&ip  sin  A<J>  + A6  cos  A4> 

-sinA<J> 

cos 

A<(> 

C = 
= 


1 A8  sin  A<f> 
0 cos  A$ 

0 sin  A<}> 


A<ji  cos  A$ 
-sin  A$ 
cos  Acf> 


(2-6) 


(2-7) 


Here  A , is  the  matrix  which  may  be  used  to  transform  the  x -,  y - and 
«=r/t  u 

z^-components  of  a vector  quantity  into  x-,  y-,  and  z-components  and  £ 


the  matrix  such  that 
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12 


y 


G X T 

Er  = [ -mg  sin  6g  0 mg  cos  <$G]  + At , F^ 

t r 

ft  = [ft  0 0]T 

r 

T = [T  0 0]T 

-s  x 
r 

Also,  in  Eq.  (2-9)  and  throughout  this  report,  a tilde  (~)  above  a 
3x1  matrix  denotes  a particular  3x3  skew-symmetric  matrix  formed  from 
the  elements  of  the  3x1  matrix.  For  example, 


0L 

rr/t 


0 -r  q 

r 0 -p 

-q  p 0 


The  equation  which  governs  the  translational  motion  of  the  rocket's 
center  of  mass  may  also  be  obtained  by  using  the  methods  of  Ref.  2.  It  is 


where 


0 

0 


h ^t/T  -j-T  £c/A  +^'" 


+ i Fex 
m -r 
t 


=t/T 

T T 


|*  + a-r/tJc/A 


(2-10) 


r„ , . = A . to  . & + A , £, 

-C/A  =t/r  -r/t  -a  =t/r  -a  . 

t r r r 


Equations  (2-9)  and  (2-10)  include  the  effects  of  the  constraint  forces 


which  cause  the  point  a to  move  with  the  turret.  They,  along  with  Eq . (2-8), 
comprise  the  equations  of  motion  for  the  tipoff  part  of  the  launch  phase. 
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2.4  Initial  Conditions  for  the  Flight  Phase 

The  main  purpose  of  the  launch  phase  equations  is  to  provide  realistic 
initial  conditions  for  the  flight  phase;  i.e.,  that  phase  beginning  when 
the  rocket  is  free  of  the  launcher's  constraints.  At  the  end  of  the  launch 

phase  simulation  in  which  the  requisite  differential  equations  are  integrated 

* 

numerically,  we  have  the  values  of  xt  = + distance  from  the  aft  end 

of  the  rocket  to  the  center  of  mass,  A<}>,  A0,  A^,  p,  q,  r and  >G  = wT(t  -tQ) , 

where  t . _ is  the  time  the  launch  phase  ends  and  t is  the  time  the  rocket 
exit  v o 

starts  to  move. 

The  initial  conditions  we  need  for  the  flight  phase  is  the  following: 
the  xG~,  y— and  zG-coordinates  of  the  rocket's  center  of  mass  C;  the  Euler 
angles  <^,  0 and  <f>  which  are  used  to  define  the  orientation  of  the  rocket- 
fixed  coordinate  system  Cxyz  with  respect  to  an  earth-fixed  system  OXYZ;  the 
inertial  components  of  the  velocity  of  C,  U,  V and  W;  P,  Q and  R,  the  x-, 
y-,  and  z-components,  respectively,  of  the  angular  velocity  of  the  Cxyz 
system  relative  to  the  OXYZ  system.  We  shall  use  a small  e as  a subscript 
here  to  denote  the  value  at  the  instant  of  exit.  Also,  the  vector  r , from 
point  G to  point  A is  introduced  at  this  point.  Using  previous  results 
and  resorting  to  Figures  2-3  and  2-4,  we  obtain  first 


(P  Q R)*  = [(p  q r)T+  Ar^T  wT]e  , 


(2-11) 


(U  V W)„  - A ou  R + k.  r, 


e £r/T  5a  + -r/t  £cM 


(2-12) 


(XG  ^ ZG)  " =t/T  (|A  ‘ + £cM 


(2-13) 


l 


Figure  2-3.  Orientation  of  the  Rocket-Fixed  System. 


To  determine  the  initial  values  of  the  Euler  angles  tp,  0 and  $>,  we 

construct  the  direction  cosine  matrix  for  the  G xG  z G system  (the  same 

matrix  as  that  for  the  Ax  y z system)  using  first  the  angles  X , 6 , Ai|», 

t t t G G 

A0,  A<J>.  This  matrix  we  call  L^j,,  where  b/I  denotes  that  the  rotation  is 
from  the  inertial  OXYZ  system  to  the  body-fixed  system  Cxyz . This  matrix 
can  also  be  constructed  by  using  the  Euler  angles  i|i,  0 and  <j>.  By  comparing 
terms  in  the  two  different  expressions  for  we  therefore  find  that 


1 


(2-14a) 
e 

(2-14b) 

<f> 

e 

sin  0g  = [sin  6G  + A0  sin  6G]g  (2-14c) 


tan  \p 


cos  <5  sin  A + Aifi  cos  A - A0  sin  6 sin  A 
hr It G G G 

cos  6n  cos  A„  - Ai[)  sin  A„  - A0  sin  6_  cos  \n 
G G G G G 


tan 


sin  6 [ Aip  cos  Aiji  - A0  sin  A(£]  + cos  <5_  sin  A if 

G G 

-sin  6_  [Aij/  sin  A4>  + A0  cos  Aijj]  + cos  cos  A 

. u b 


We  have,  at  this  point,  developed  the  basis  for  simulating  the  launch 
phase  and  have  also  obtained  equations  from  which  we  can  determine  initial 
conditions  for  the  flight  phase,  given  the  final  conditions  for  the  launch 
phase.  In  the  next  section, we  develope  the  flight  phase  model. 


SECTION  3.  FLIGHT  DYNAMICS 


\ 


\ 
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3.1  Introductory  Comments 

In  this  section,  we  develop  the  mathematical  model  for  the  rocket. 

This  model  is  essentially  the  same  as  that  developed  in  Ref.  2 and  used 
there  as  the  basis  for  the  simulation  of  the  free-flight  phase  of  a "free," 
or  uncontrolled,  rocket.  The  differences  in  the  two  models  are  that 
(1)  additional  forces  and  moments  due  to  the  deflection  of  the  sustainer 
motor's  thrust  vector  are  included  here  and  (2)  the  time  variations  of  the 
mass  of  the  rocket  and  of  its  centroidal  moments  of  inertia  are  assumed, 
for  the  purposes  of  this  study  to  be  prescribed  a_  priori  as  explicit  func- 
tions of  time. 

3.2  Physical  Model  for  the  Rocket 

The  rocket  is  modeled  physically  as  a variable  mass  body  composed  of 
a rigid  portion  and  a fluid  portion  (see  Fig.  3-1).  The  rigid  portion 
consists  of  the  unexpendable  parts  of  the  rocket,  plus  the  part  of  the 
solid  propellant  not  yet  expended  at  time  t.  The  fluid  portion  consists 
of  that  part  of  the  solid  propellant  which  has  been  burned  (or  is  burning) 
and  is  flowing  as  a fluid  within  the  rocket,  having  not  yet  passed  through 
the  rocket  nozzles.  For  the  purposes  of  modeling  the  internal  flow,  the 
geometry  of  the  rocket  is  simplified  by  assuming  that  the  flow  is  within 
the  core  of  a cylindrical,  solid  propellant  charge.  Since  the  internal  flow 
model  is  actually  used  only  in  estimating  the  jet-damping  torque,  this  simple 
model  is  considered  sufficient. 

The  rocket,  apart  from  the  solid  propellant,  consists  (see  Fig.  3-2) 
first  of  all  a shell,  or  fuselage,  which  contains  the  warhead,  the  solid 
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Solid,  Non-expendable  Part 


Fluid  Part 


Y//////////////////SS// VY//////////////,y////////YY/////////Y//. 


^Solid  Expendable  Part 


Figure  3-1.  Solid  and  Fluid  Portions  of  the  Rocket. 
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propellant,  the  rocket  nozzles,  a pitch  rate  gyro,  a roll  angle  gyro,  the 
actuating  mechanism  for  the  thrust  deflector  vanes,  the  vanes  themselves, 
the  electronic  equipment  needed  to  receive  and  process  guidance  commands 
and  a device  which  controls  the  motion  of  destabilizing  fins.  The  last 
six  items  listed  will  be  discussed  more  fully  in  the  section  on  guidance 
and  control  (see  Section  4).  Four  wings  are  attached  to  the  rocket's 
fuselage  in  a cruciform  configuration.  These  wings  fold  out  as  the  rocket 
exits  the  launch  tube.  Also,  there  are  four  small  fins  which  are  extended 
through  the  slots  near  the  nose  of  the  rocket.  They  are  the  destabilizing 
fins  mentioned  above.  They  are  extended  in  a pre-programmed  manner  so  that 
the  aerodynamic  center  of  the  rocket  moves  forward  as  the  center  of  mass 
moves  forward  (due  to  propellant  consumption),  with  the  result  that  a fairly 
constant  static  margin  is  maintained  during  the  rocket's  powered  flight. 

3.3  Equations  of  Motion 

The  equations  of  motion  for  the  rocket  may  be  obtained  using  the  pro- 
cedure followed  in  Ref.  1,  or,  for  example,  that  of  Ref.  3. 

Translational  Motion 

The  equation  which  governs  the  translation  of  the  center  of  mass  of 
the  rocket  may  be  written  in  the  form. 


i V,  = F.  + F + F_  , + F „ , + F , 

-C  -T  -c  -Coriolis  -unsteady  -ex’ 


(3-1) 


where  F is  the  thrust  "force"  (see  Ref.  4)  which  is  generated  by  the  outflow 


of  the  fluid  part  of  the  rocket,  F is  the  control  force  (also  generated 

—c 


by  the  same  mechanism  as  F , but  segregated  here  for  convenience).  Ecoriolis 


The  word  "force"  is  enclosed  in  quotation  marks  here  to  indicate  that  it 
pertains  to  a term  in  the  equation  of  motion  which  is  treated  as  a force,  but 
actually  arises  due  to  the  fact  that  mass  is  being  expelled  from  the  rocket 

and  strictly  speaking  is  a reaction  force  on  the  solid  part  of  the  rocket  due 


hJ 
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is  the  "force"  required  to  accelerate  the  mass  of  fluid  within  the  rocket 

transversely,  F , is  the  "force"  resulting  from  the  time  rate  of 

-unsteady 

change  of  fluid  flow  properties  at  a given  spatial  position  and  Fgx 
includes  all  forces  on  the  rocket  due  to  external  influences  such  as  the 
atmosphere  and  gravity. 

In  most  cases,  F,,  . , . and  F . , are  not  significant.  Hence, 

we  will  neglect  these  two  "forces."  Then,  by  writing 

V=Ui+Vj+Wk,  (3-2) 

~C 

A A A 

where  i,  j and  k are  unit  vectors  fixed  to  the  rolling  rocket  which  has 
angular  velocity, 

2=Pi+Qj+Rk,  (3-3) 

O • ~ A • a 

and  letting  yc=  Ui+Vj+Wk,  we  may  rewrite  Eq.  (3-1)  in  the  form, 

V = - x V + - [F  + F + F ] (3-4) 

-C  - - m -T  -c  -ex 

If  the  rocket  is  imperfect,  which  is  usually  the  case,  then  to  account  for 
thrust  vector  misalignment  we  write 

F_  = T cos  a cos  a i + T cos  a sin  a j - T sin  a k,  (3-5) 
-T  yz  y z J y ’ 

where  T is  the  magnitude  of  the  thrust  and  and  ot^  are  constant  angles 

(as  shown  in  Fig.  3-3  ) which  are  measures  of  the  amount  and  direction  of 

the  mechanical  thrust  misalignment  present.  It  should  be  pointed  out  that 

these  angles  represent  the  "equivalent  thrust  misalignment"  resulting  from 

the  combined  mechanical  misalignment  of  both  the  booster  motor  nozzles. 


y 
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The  control  force  is  discussed  more  fully  in  Section  4,  but  here  we 
note  that  (see  Fig.  3-4). 


Fc  = +F[-sin  6 j + cos  6 k]. 


(3-6) 


where  F is  the  magnitude  of  the  control  force  and  6 is  a constant  angle 
representing  the  rotation  of  the  control  plane  with  respect  to  the  xz- 
plane . 

The  external  force  which  acts  on  the  rocket  is  the  resultant  of  the 
aerodynamic  force  and  the  force  due  to  gravity,  hence. 


where 


and 


F = F.  + F , 
-ex  -A  -g  ’ 


FA  = -A  i - Y j - N k 


Fg  - mg[-sin  6 i + cos  6 sin  <j>  j + cos  6 cos  (j>  k]  . 


(3-7) 


(3-8) 


(3-9) 


In  Eq.  (3-8),  A is  the  axial  aerodynamic  force  component,  Y is  the  y- 
component  of  the  aerodynamic  force  and  N is  the  normal  aerodynamic  force 
component  (see  Fig.  3-5).  In  Eq.  (3-9),  0 and  £ are  two  of  the  Euler 
angles  which  are  used  to  define  the  orientation  of  the  rocket-fixed  coordinate 
system  Cxyz  (see  Fig.  2-3). 

For  the  purposes  of  this  study,  the  aerodynamic  forces  are  assumed  to 
be  expressible  as 


Control  Axis 


Figure  3-4,  Control  Force  Orientation. 


Normal  Force 

-Nfc 


Yawing  Moment 


Figure  3-5,  Aerodynamic  Forces  and  Moment 


where  p is  the  atmospheric  density,  S is  the  reference  area  upon  which 

the  aerodynamic  coefficients  C.  , C and  C„  are  based  and  V„  is  the  absolute 

A y N R 

value  of  the  relative  wind  that  the  rocket  encounters.  Herein 


VR  =V(UW-U)2  + (Vw-V)2  + (Ww-W): 


(3-10) 


where  U^,  and  are  the  x-,  y-  and  z-components  of  the  velocity  of  the 
wind . 

The  usual  assumption  that  is  a function  of  angle  of  incidence  and 
Mach  number,  M,  is  made  and  it  is  further  assumed  that  the  incidence  angle + 
is  small.  Then,  the  following  forms  for  C^,  and  CN  are  adopted: 


CA  " 


(3-lla) 


Cy  = -S(M)  3 

a 


(3-llb) 


a » 


(3-llc) 


where  C (M)  and  CM(M)  are  tabulated  functions  and  a = tan-1[(W-W  )/(U-U  )] 

A N W W 

o a 

and  3 = sin-1  ^(V-Vw)/VR]  are  t*ie  aeroc*ynamic  angles  of  attack  and  sideslip, 

respectively. 

In  the  simulation  program,  the  differential  equations  for  V and  W are 
not  integrated.  Instead,  the  variables  = tan-1 (W/U)  and  Bx  = sin_1(V/Vc) 
are  used  in  place  of  V and  W for  computational  purposes.  The  variables  a^and 

Bxare  the  aerodynamic  angles  of  attack  and  sideslip,  respectively  if_  there 

• • • 

is  no  wind.  The  derivatives  V and  K are  computed  using  Eq.  (3-4)  and  a^and 

B^are  obtained  from 


a = (UW  - WU)/(U2  + W2) 


(3-12a) 


e = [V(U2  + W2)  - V(UU  + WW)]/[V^Aj2+W2  ] 


(3-12b) 


The  kinematic  equation  for  translation  is  obtained  by  noting  that  the 
velocity  of  the  center  of  mass  of  the  rocket  may  be  written  either  as 


Vc=  XI+YJ  + ZK, 


(3-13) 


where  I,  J and  K are  unit  vectors  associated  with  the  earth-fixed  OXYZ  system, 
or  as 


V = U i + V j W k , 
-C 


(3-14) 


Since  the  components  of  in  the  two  systems  are  related  by  the  transforma- 
tion matrix  (see  Fig.  2-3)  used  previously  in  Section  2, 


cdc<p 

c6s<p 

-s8 

-c<psip 

c4>c^ 

s q>c  0 

+ s<psdc\p 

+ S<J)S0SlJj 

s<ps>p 

— S<J>C^ 

cd>c0 

+ c<j>s©ct// 

+ Ctpsdstp 

(3-15) 


The  required  kinematic  equation  is 


=l/b  V ’ 


(3-16) 


where  LT  ,,  = L,  /T 
=*I/b  =b/I 
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Rotational  Motion 


The  equation  governing  the  rotational  motion  of  the  rocket  about  its 
center  of  mass  C is 


I il  + 8 X I • a • T.  + T + T ^ , + T„  ,nj  + T 
= - - = - -T  -c  -unsteady  -Coriolis  -ex 


(3-17) 


where  I is  the  centroidal  inertia  dyadic  of  the  rocket,  is  the  "torque"^ 

about  C due  to  the  thrust,  T is  the  control  "torque"  about  C,  T . is 

-c  -unsteady 

the  "torque"  about  C due  to  the  unsteadiness  of  the  flow,  T_  . , . is  the 
"torque"  about  which  is  due  to  the  time  rate  of  change  of  the  direction  of 
the  motion  of  the  fluid  which  is  flowing  inside  the  rocket  caused  by  rotation 
of  the  rocket  and  is  the  torque  about  C due  to  aerodynamic  reactions. 

The  inertia  dyadic  I is  assumed  to  be  general  in  form  to  account  for 
mass  imbalance  which  may  be  present  due  to  manufacturing  imperfections.  The 
matrix  counterpart  of  I expressed  in  terms  of  rocket-fixed  elements  is 


I 

-I 

-I 

XX 

xy 

xz 

I = 

-I 

I 

-1 

(3-18) 
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r 

xy 

yy 

yz 

-I 

-i 

I 

xz 

yz 

zz 

the  ideal  case.  I 

= I 

I 

0 and  I - I 

xy 

xz 

yz 

yy  zz 

The  torque 

due  to 

thrust 

is  expressed  explicitly  as 

ii 

H 

HI 

<*T  ' H> 

x Ft  , 

(3-19) 

+The  word  "torque"  is  enclosed  in  quotation  marks  here  to  indicate  that 
it  pertains  to  a term  in  the  equation  of  motion  which  is  treated  as  an  external 
torque,  but  actually  arises  from  the  fact  that  mass  is  being  expelled  from 
the  rocket  and  strictly  speaking  is  a reaction  torque  on  the  solid  part  of 
the  rocket  due  to  the  internal  flow  or  expelled  mass.  The  quotation  marks 
will  be  dropped  in  later  references  to  the  same  term. 


I 

1 


1 

j 


1 
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where  S.^,  is  a vector  from  an  arbitrary  point  in  the  non-expendable  part 
of  the  rocket  to  a point  on  the  line  of  action  of  the  resultant  thrust 
vector.  Also,  in  Eq.  (3-18)  u is  a vector  from  a to  C (see  Fig.  3-6). 

The  control  torque,  like  the  control  force,  will  be  discussed  more 
fully  in  Section  4.  However,  we  note  here  that  it  has  the  form. 


Ic  -<il  - 2>  * lc  • 


0-20) 


where  £T  is  a vector  from  point  a to  the  line  of  action  of  the  control 
force. 

For  the  purposes  of  this  study  the  torque  due  to  the  unsteadiness  of 
the  flow  is  assumed  to  be  negligible  compared  to  other  torques. 

The  Coriolis  torque  is  the  so-called  "jet-damping"  torque,  TjD>  and 
is  retained.  By  assuming  a cylindrical  solid  propellant  charge  and  constant 
mass  flow  rate  per  unit  length  of  the  change,  an  approximate  expression  for 
^Coriolis  may  °btainec*-2  The  expression  used  in  this  simulation  study  is 


T = T„  . ..  = — 2 CTTJ0  j + R k] 

-JD  -Coriolis  JD  ^ J 1 


(3-21) 


Here,  the  jet  damping  coefficient,  CjD>  is  given  by 

CJD  = - m[(ux  - x'/3)*p  +(xp2/3]  , (3-22) 

where  m is  the  time  rate  of  change  of  the  rocket's  mass,  u^  = u • i,  £ 
is  the  length  of  the  propellant  charge  and  xj^  and  xJ,  are  the  coordinates 
of  the  aft  and  fore  ends,  respectively,  of  the  propellant  charge  as  measured 
from  point  a as  shown  in  Fig.  3-7. 
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The  external  torque  about  the  rocket's  center  of  mass  Is  due  to  the 
resultant  aerodynamic  force.  This  aerodynamic  torque  may  be  expressed  in 


the  conventional  form, 


T =T=Li  + Mj+Nk 
-ex  -A  J 


(3-23) 


where  L,  M and  N are  the  rolling,  pitching  and  yawing  components,  respec- 
tively, of  T . The  components  L,  M and  N may  be  expressed  in  the  forms, 


L = 1/2  p s vj  d Cj 

M = 1/2  p S d C 
R m 

N = 1/2  p S V*  d C , 
R n 


(3-24a) 

(324b) 


(3-24c) 


where  C£,  and  are  the  usual  aerodynamic  coefficients  and  d is  the 
reference  length  upon  which  these  coefficients  are  based. 

In  general. 


c t = C£(M,o,e,P) 

C = C (M,a,a,Q) 
m m ’ * ’ 


(3-25a) 

(3-25b) 


Cn  = cn(M’B’6’R)’ 


(3-25c) 


where  M is  used  here  to  denote  the  flight  Mach  number.  For  small  a ^ and  6^, 
which  we  assume  is  the  case,  can  be  adequately  represented  by  the  ex- 
pression. 


Cz  = [d/(2VR)]C£  (M)  P, 


(3-26a) 


Note  that  the  symbols  M and  N are  used  elsewhere  for  the  flight  Mach 
number  and  the  normal  component  of  the  aerodynamic  force,  respectively. 

This  notation  is  also  conventional  and  it  should  be  clear  from  *"he  context  in 
which  the  symbols  appear  which  meaning  is  to  be  attached  to  each  of  them. 
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where  is  a non-dimensional  stability  derivative  defined  as 


c*  E (2VR/d)(9q/3P)p=0 

P 


Also,  for  small  a,  8,  P and  Q, 


where 


C = C a + [d/(2V  )]  C Q 

mm  K m 7 

a q 


c = (3C  /3a) 
m m _ 

a a=0 


C = (3C  /3Q)(2VB/d) 

“q  m R Q=(/ 


(3-26b) 


Since  the  rocket  is  geometrically  symmetiic,  C = - C and  C = C so 

n„  m n m , 

8 a r q’ 

that  we  also  write 


C * - C 6 + [d/(2V  )]C  R . 
n m R m 

a q 


(3-25c) 


In  this  report,  the  quantities  C.  , C and  C are  assumed  to  be  given 

Si  m m e 

pa  q 

functions  of  the  flight  Mach  number  in  the  form  of  tabulated  values.  The 

stability  derivative  C is  a function  of  the  static  margin;  i.e., 

a 

(x  - x )/d  where  x and  x are  the  distances  from  the  rocket's  nose  to 
cp  eg  cp  eg 

its  aerodynamic  center  and  its  center  of  mass,  respectively.  This  quantity 


is  a function  of  M because  the  distance  x varies  directly  with  Mach  number 

cp 

and  x£g  varies "indirectly  with  Mach  number 7 since  as  the  rocket's  fuel  is 

consumed  and  its  center  of  mass  moves  forward,  the  Mach  number  increases  until 

burnout  of  the  booster  motor.  The  destabilizing  fins  shown  in  Fig.  3-2  are 

extended  as  the  fuel  is  consumed  in  such  a way  that(x  - x ) remains  fairly 

cp  eg 

constant  during  powered  flight;  hence,  C also  remains  fairly  constant.  The 

a 
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destabilizing  fins  do_  not  make  the  rocket  unstable,  but  they  only  make  it 
t less  statically  stable  than  it  would  otherwise  be  at  Mach  numbers  near  and 

above  unity. 

The  basic  mathematical  model  of  the  rocket  has  been  presented.  We 
turn,  in  the  next  section,  to  the  task  of  formulating  guidance  and  control 
algorithm  for  the  system  composed  of  the  rocket  and  the  sensing  and  guidance 
command  generating  equipment  on  the  launcher. 
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SECTION  4.  GUIDANCE  AND  CONTROL 


4.1  General  Comments 

The  guidance  and  control  of  the  rocket  prior  to  radar  acquisition  are 
considered  in  this  section.  First,  the  problem  of  guiding  the  rocket  toward 
a moving  target  is  addressed.  Then,  the  control  algorithm  used  in  this  in- 
vestigation is  described. 

4 . 2 Guidance 

It  is  assumed  that  guidance  commands  are  formulated  in  the  launch 
vehicle  by  means  of  a special  purpose  guidance  computer.  The  input  to  the 
guidance  computer  is  angular  data  acquired  by  the  launcher-fixed  IR  sensor. 
This  information  is  available,  however,  only  if  the  aft  end  of  the  rocket 
(which  contains  flares)  is  within  the  Jield  of  view  (FOV)  of  the  IR  sensor. 
The  sensor  is  assumed  to  have  two  FOV's  - one  for  acquisition  of  the  rocket 
and  a more  narrow  one  for  use  after  acquisition  by  the  sensor,  but  prior  to 
radar  acquisition  (see  Fig.  4-1),  The  guidance  computer  is  assumed  to  con- 
tain a pre-programmed  range  versus  time  relation,  so  that  the  spherical 
coordinates,  (R  , e , e ) (see  Fig.  4-1)  where  R„  is  the  range  of  the  rocket 
and  z^  and  z^  are  the  "aximuth"  and  "elevation"  angles,  respectively,  of  the 
rocket  with  respect  to  Gx  y z coordinate  frame.  Furthermore,  for  the 

0 u u 

purposes  of  this  study,  it  is  assumed  that  values  of  R„,  e and  e are 

G z y 

available.  These  may  be  estimated  by  using  the  range  versus  time  relation- 
ship, the  angular  measurements  and  the  gimbal  rates  of  the  radar  (which  is 


assumed  to  perfectly  track  the  target) . For  use  in  the  current  control 
algorithm,  the  linear  devi.-.tions  of  the  rocket  from  the  line  of  sight  (LOS) 
from  the  launcher  radar  to  the  target  (This  LOS  is  the  extended  x^-axis) 
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are  used.  These  deviations  are  and  z^  which  may  be  computed  by  using 


the  equations. 


y„  = R.  cos  e sin  e 
G G y z 


(4-la) 


zG  = - rg  sin  ey 


(4-lb) 


The  time  rates  of  change  of  y_  and  z„  are  also  used  in  the  control  algorithm. 

Cj  v 

These  may  be  obtained  from 


y_  = R„  cos  e sin  e - R_  e sin  e sin 
G G y zGy  y z 


+ R„  e cos  e cos  e 
G z y z 


(4-2a) 


z = — [R_  sin  e + R„  e cos  e ] 
G G y G y y 


(4-2b) 


In  an  actual  system,  these  rates  could  be  estimated  by  using 


yG  = ' RG  XG  C0S  6G  + RG  £z 


(4-3a) 


ZG  = - RG  6G  ' RG  ey 


(4-3b) 


for  small  ey,  Ey,  Eg  and  pre-radar  acquisition  values  of  Rg. 

In  this  report,  we  consider  only  constant-speed  targets  (see  Fig.  4-2). 
For  such  targets,  the  angles  x„  and  6 are  given  by 

O Kj 


>c  * tan“1(yT/xT) 


(4-4a) 


<5g  = sin-1  (-zT/pT) 


(4-4b) 


respectively,  where 


T 

ST 


/Xi 


+ x (t-t  ) , 

(4-5a) 

T o 

+ yT(t-to)’ 

(4-5b) 

+ zT(t-to), 

(4-5c) 

+ ?T  + ZT 

(4-5d) 

and  x,j, , yT  and  z T are  the  coordinates  of  the  target  in  a non-rotating 
coordinate  system  with  origin  at  G.  It  follows  from  Eqs.  (4-4)  that 


~ (y-f/Xj,  - yT  xt/xT)cos2  XG  (4-6a) 

and 

- C-zt/pt  + zT  p£/p|)/cos  «G  (4 -6b) 

where  pt  = (XT  + yT  yT  + z T zt)/pt- 
4.3  Control 

Control  of  the  rocket  is  physically  accomplished,  as  mentioned 
previously,  by  causing  the  deflector  vanes  to  dwell  in  either  a positive  or 
negative  position,  thereby  producing  a torque  about  the  control  axis 
(see  Fig.  3-4)  of  the  rocket.  This  torque  produces  an  angular  acceleration 
about  either  the  positive  or  negative  control  axis.  The  integrated  effects 
of  the  torque  are,  first  of  all,  an  angular  rate  and,  secondly,  an  attitude 
change.  The  latter  results  in  a change  of  the  rocket’s  angle  of  incidence 
and  an  accompanying  change  in  the  aerodynamic  force  and  torque  on  the  rocket. 
Finally,  the  altered  aerodynamic  force  produces  a change  in  the  rocket's 
flight  path.  The  duration  of  the  dwell  and  the  roll  orientation  of  the 
rocket  while  it  is  taking  place  determine  the  resultant  change  in  its  flight 
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The  purpose  of  the  control  system  is  two-fold.  The  first  part  of  its 
purpose  is  to  "stabili2e"  the  rocket’s  angular  velocity  in  such  a way  that 
angular  rates  due  to  disturbances,  and  not  guidance  commands,  are  nulled 
out  rapidly.  Secondly,  it  must  produce  the  attitude  changes  dictated  by  the 
guidance  commands. 

Since  the  control  force  is  either  "on"  or  "off",  the  control  system  is 
of  the  bang-bang  type.  Thus,  in  formulating  a control  algorithm  the  questions 
of  when  to  turn  the  control  force  on  (and  in  what  direction)  and  how  long  to 
leave  it  on  (the  dwell  time)  must  be  answered. 

Control  Algorithm 

We  let  | At  j and  j At  1 denote  the  dwell  times  of  the  deflector  vanes 

y 2 

for  control  of  the  rocket’s  motion  in  the  XgzG~  and  xGyG~planes,  respectively. 

The  signs  At  and  At  determine  whether  F = + F , as  will  be  explained 
y C Cmax 

in  what  follows.  We  also  let  5^  and  6z  denote  the  angles  through  which  the 

rocket  rolls  during  a "longitudinal"  (x_z_,-plane) , or  "lateral"  (x„y_-plane) , 

control  period.  To  be  more  specific,  we  assume  that  positive  control 

torques  are  desired.  Then,  to  control  longitudinal  motion,  the  control 

force  is  given  by  (see  Fig.  4-3) 

F * F , if  2mr  - 6 (1  < 6 + 6 < 2nn  + 6 /2  , (4-7a) 

c c ’ y * y 

max 

and  by 

F - - F , if  (2n+l)7r  - 6 /2  < $ + 6 (2n+l)ir  + 6/2.  (4-7b) 

cmax  y y 

Here  n * 0,1,2,....  Similarly,  for  lateral  control, 

F = F , if  [(4n+l)/2]u  - 6^/2  £ $ + 6 <_  [(4n+l)/2]n  + 6z/2  , 
cmax 


(4-8a) 


1 
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and 

Fc  = - Fc  , if  [(4n+3)/2]ir  - 6/2  < <j>  + 6 <_  [ (4n+3) /2 ]tt  + 6/2. 
cmax  z 

(4-8b) 

Ideally,  tne  control  forces  given  by  Eqs.  (4-7)  and  (4-8)  act  as  space- 

fixed  forces  which  produce  space-fixed  torques  and  the  latter  result  in 

changes  in  the  rocket’s  attitude.  Since  the  rocket  is  rolling,  there  is 

some  coupling  between  its  pitch  and  yaw  motions.  But,  since  it  is  also 

"slender"  (i.e.,  1^  <<  I - Izz)»  this  coupling  is  not  severe. 

The  roll  angle  $ must  be  known  if  is  to  be  generated  properly.  This 

angle  is  assumed  to  be  measured  by  a roll  gyro  on  board  the  rocket . 

The  dwell  times  At  and  At  are  determined  from  the  control  laws 
y Z 

4ty  * ' K1  - V > <4-9a> 


and 


At 


" " Ki  (rnr  ” RNR0) ' 


(4-9b) 


where  1C  is  a constant,  Q is  the  angular  rate  of  the  rocket  about  the  y — 

axis  [which  is  one  of  the  three  axes  of  the  non-rolling  Cx^y^z^  coordinate 

system  (see  Fig.  4-4)],  R^R  is  the  angular  rate  of  the  rocket  about  the 

Xj^-axis  and  and  R^R  are  commanded  angular  rates  about  the  y - and 

c c 

z^-axes,  respectively.  The  angular  rates  and  R^  are  computed  from 


Qho  * Q cos  d>  — R sin 

NR 


(4-10a) 


and 


RjjR  = Q sin  $ + R cos  #, 


(4-10b) 


respectively,  Also,  the  commanded  angular  rates  are  determined  using  the 


equations 


QNR  k.  ZG  + k,  ZG 
cl  1 


_ _ ^2  ^3  • 

Sr  " k.  yG  k,  yG 
c 1 1 


(4-llb) 


respectively,  where  and  are  time  varying  gains  which  have  the  forms. 


K2  K20  to/t 


(4-12a) 


K30  ‘ K30  '</«  , (t  i ‘o'  • 


(4-12b) 


Here,  Kjq  and  K^q  are  constant  and  tQ  is  the  time  the  control  system  is 
enabled  (This  is  a short  time  after  the  rocket  exits  the  tube  and  corresponds 
to  the  time  the  sustainer  motor  is  ignited.) 

We  note  that,  for  an  actual  single-axis  system,  both  Q and  R may  not 
be  measured  by  rate  gyros.  In  fact,  there  will  probably  be  only  one  rate 
gyro  which  measures  only  the  angular  rate  of  the  rocket  about  the  control 
axis.  However,  6y  and  6^  are  herein  restricted  to  small  angles,  so  that  the 
value  of  (R^)  used  to  calculate  At  (Atz)  is  approximately  the  angular 
rate  about  the  control  axis  just  prior  to  a control  phase. 

In  an  actual  system,  values  of  and  R^  are  calculated  by  the 

guidance  computer  and  sent  to  the  rocket  via  a radio  link.  These  values  are 
used  by  the  rocket's  on-board  systems  to  generate  At  and  At  commands.  In 

y z 

the  present  simulation,  the  procedure  used  is:  (1)  Specify  priori 

6 - 6 and  an  associated  At  =6  /£2,  where  £5  is  the  average  roll 

y z max  y 

max  max  max 

rate  of  the  rocket.  (2)  Compute  Aty  (Atz)  and  check  to  see  if  it  is  less 

than  At  . (3)  If  At  < At  (At  < At  ),  then  compute  6 “ | At  |ft  (6  ■ 

max  y — max  z — max  y ' y'  z 

|AtzfG).  (4)  If  Aty  > AtMx  (Atz  > Atmx),  then  set  Aty  = sign  (At^At^ 


(At  * sign  (At  )At  ) before  computing  5 (6  ) . (5)  Check  to  see  if  the 
z z m&x  y z 
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rocket's  roll  attitude  satisfies  Eqs.  (4-7)  (Eqs.  (4-8)).  (6)  If  it  does, 

control  is  exercised.  If  not,  no  action  is  taken.  (7)  If  the  time  since  the 

beginning  of  this  control  phase  is  less  than  t,  the  control  force  magnitude 

is  adjusted  so  that  a ramp  up  from  zero  to  F is  obtained  in  t seconds 

max 

during  the  beginning  of  the  phase.  Also,  if  the  time  before  the  end  of 
i this  control  phase  is  less  than,  or  equal  to,  t,  a ramp  down  in  F is 

? V- 

max 

generated  during  the  last  x seconds  of  the  phase.  (8)  If  |At^|  < 2x 
( | At z | < 2x),  no  control  is  exercised. 

Control  Gain  Selection 

An  approximate  linear  stability  analysis  was  conducted,5  assuming  that 

the  rocket's  roll  rate  is  rapid  enough  that  the  average  values  of  the  control 

torques  generated  during  one  period  of  the  rolling  motion  could  be  considered 

continuous  torques  with  respect  to  the  non-rolling  system  Cx^y^RZ^.  This 

analysis  provided  some  indication  of  the  values  of  K^,  and  K^q  which 

should  be  used, but  the  responses  predicted  by  the  linear  analysis  were  not 

the  same  as  those  obtained  via  nonlinear  simulations.  This  might  be  expected, 

since  the  control  problem  is  actually  very  nonlinear.  Hence,  preliminary 

values  of  the  gains  were  chosen  using  physical  and  hueristic  reasoning  and 

"refined"  values  obtained  by  changing  the  gains  selectively  and  generating 

trajectories  with  the  simulation  code. 

The  reasoning  used  to  pick  the  preliminary  gains  was  as  follows.  The 

commanded  transverse  angular  rates  of  the  rocket  are  assumed  to  be  limited 

by  Q . Hence,  in  one  control  period  the  control  torque  total  impulse 
Cmax 

must  be  no  more  than  that  necessary  to  produce  a rate  Qc  . For  a given 
[ max 

maximum  control  region  width,  6-6  , the  duration  of  control  application 

i ^max  zmax 

i 

L> 


1 

•1 

j 


is  small,  so  that  the  effective 


is  At  =6  /fl.  Hence,  if  6 
max  y y 

max  max 

magnitude  of  the  control  force  is  approximately  F , if  I^,  is  the  value 

max  © 

of  the  lesser  of  I or  I at  the  beginning  of  controlled  flight  and  if  l 


yy 


zz 


is  the  effective  moment  arm  of  about  C,  then  the  change  in  magnitude  of 
the  angular  rate  is  approximately  (longitudinal  control) 


Q = 1 F At  /IT 
c c c max  T 

max  max  o 


(4-13) 


Now,  from  Eq.  (4 -9a) (with  Q^R  = 0 and  = Q£  ) 

C max 


At  = 1C-  Q 
max  1 c 


max 


Hence,  from  Eqs.  (4-13)  and  (4-14),  we  have 


K1  “ h Fc  > 

o max 


(4-14) 


(4-15) 


By  replacing  z and  z with  their  expected  maximum  values  at  the  beginning 

U (j 

of  the  controlled  flight  and  by  Q , in  Eq.  (4-lla),  we  obtain 

c max 


% K20'ZG'  K30'zg' 

max  max  max 


(4-16) 


By  considering  the  results  of  the  linear  stability  analysis5  and  some  trial 
simulation  results,  it  was  found  that  equal  numerical  values  for  K^q  and  K^q 
resulted  in  reasonable  response  characteristics.  Hence,  (except  for  units) 


*20=  K30=  Qc  ^ZG 


(4-17) 


max  max 

were  chosen  as  preliminary  values  of  the  "guidance"  control  gains. 

As  stated  above,  the  preliminary  values  of  K^,  K^q  and  K^q  were 
adjusted  to  obtain  values  which  result  in  reasonable  response  characteris- 


tics of  the  typical  rocket  simulated  in  this  study. 
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5.  SIMULATION  RESULTS 

5.1  General  Comments 

To  check  out  the  operation  of  the  simulation  codes  and  to  study  the 
effects  of  the  various  anomalies,  such  as  winds,  several  simulations  were 
carried  out.  In  the  course  of  these  simulations,  several  different  numerical 
integration  algorithms  were  used  to  solve  the  flight  phase  equations  of 
motion,  while  the  launch  phase  equations  were  integrated  using  only  a 
standard  fourth-order  Runge-Kutta  algorithm.  The  procedures  used  to  solve 
the  flight  phase  equations  were  (1)  a variable-step,  fourth-order,  Runge- 
Kutta  algorithm,6  (2)  a variable-step  Hamming  predictor-corrector  algorithm,6 
(3)  a fixed-step,  fourth-order,  Runge-Kutta  algorithm  and  (4)  a fixed-step 
Euler  algorithm. 

The  variable-step  algorithms  were  not  satisfactory  because,  during  time 
intervals  in  which  the  control  force  changes  very  rapidly,  the  step-size 
adjusting  segments  of  the  algorithms  failed  to  find  a suitable  step-size. 

This  occurred,  not  because  of  inaccuracies  in  the  calculations,  but  as  a 
result  of  the  methods  used  to  check  on  the  accuracy  of  the  computations. 

The  fixed-step  algorithms  were  tested  by  using  various  step  sizes  and 
a "best"  step  size  was  obtained  for  each.  These  "best"  step  sizes  are  de- 
fined as  those  which  produced  consistent  solutions  and  which  are  such  that 
step  sizes  larger  or  smaller  produced  inconsistent  results. 

Simulations  using  the  fixed-step,  Runge-Kutta  algorithm  with  its  "best" 
step  size  require  more  computer  time  than  their  Euler  algorithm  counter- 
parts. Hence,  the  fixed-step,  Euler  algorithm  was  used  to  generate  the 
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results  presented  herein.  However,  the  simulation  code  contains  both 
algorithms  and  the  one  used  is  optional  with  the  user. 

5.2  Data  Used  in  the  Simulations 
Most  of  the  data  for  the  launch  phase  simulations  is  given  in  Table 

5-1.  The  remaining  data  is  the  angular  rate  of  the  launcher  turret  and 
the  tipoff  distance.  These  were  different  for  each  of  the  cases  for  which 
results  are  given  and  may  be  found  in  Table  5-2. 

The  data  for  the  flight  phase  simulations  is  given  in  Table  5-3  and 
Appendix  C.  The  physical  characteristics  of  the  rocket  are  given  in 
Table  5-3.  The  manner  in  which  aerodynamic  parameters  were  calculated 
is  explained  in  Appendix  C.  Values  used  for  the  control  gains  were 
«*  0.125  sec2,  = 0.03  rad-sec/m  and  K^q  = 0.03  rad-sec2/m. 

5.3  Simulation  Results 

Since  the  question  of  whether  the  rocket  is  acquired  by  the  IR  sensor 
and  guided  in  such  a manner  that  radar  acquisition  is  possible  is  the  one 
of  interest,  only  the  flight  phase  results  of  the  simulations  are  presented 
here.  The  results  of  the  corresponding  launch  phase  simulation  provided  the 
initial  conditions  for  each  flight  phase  simulation.  The  results  for  ten 
simulations  are  given  here.  These  may  be  classified  as  follows: 

Simulation  1 - No  Anomalies 
Simulation  2 - Tipoff 
Simulation  3 - Mass  Imbalance 
Simulation  4 - Thrust  Misalignment 
Simulation  5 - Cross  Wind 

Simulations  6 and  7 - "Positive"  Turret  Rates 
Simulations  8 and  9 - "Negative"  Turret  Rates 
Simulation  10  - "Combination"  of  Anomalies 


TABLE  5-1.  LAUNCH  PHASE  DATA 


Rocket 

LENGTH  2.4  (m) 

DIAMETER  0.16  (m) 

MASS  62.4  (kg) 


Ix  =0.325  (kg-m2) 
ly  = Tz  = 21.3  (kg-m2) 

THRUST  16850.0  Nt) 

SPIN  TORQUE  100.0  (Nt-m) 

Other  Data 

|a  = (-0.95  1.03  0.0)T  (m) 

|g  = (-0.95  0.0  -1.015)1  (m) 

5c/A(INITIAL)  = (1.001  0.0  0.0)T  (m) 

*E  - 2.4  (m) 
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TABLE  5-3.  FLIGHT  PHASE  DATA. 

REFERENCE  AREA  FOR  AERODYNAMICS  - 0.021  (m2) 

REFERENCE  LENGTH  FOR  AERODYNAMICS  - 0.16  (m) 

TIME  OF  BOOSTER  BURNOUT  - 1.65  (sec) 

MASS  - (INITIAL)  - 63.2  (kg) 

MASS  - (AT  BOOSTER  BURNOUT)  - 58.3  (kg) 

MOMENTS  OF  INERTIA  (INITIAL)  - Ix  = 0.325  (kg-m2) 

V “ I*  " 21.5  (kg-m2) 

MOMENTS  OF  INERTIA  ( AT  BOOSTER  BURNOUT)  - Ix  - 0.299  (kg-m2) 

I -I  =15.46  (kg-m2) 

y z 

u (INITIAL)  - (1.001  0.0  0.0)T  (m) 

U (AT  BOOSTER  BURNOUT)  - (1.110  0.0  0.0) T (m) 

£t  - (0.124  0.0  0,0)T  (m) 

|x*|-  0.158  (m) 

- 0.793  (m) 

SLAT  - 0,0  (m) 

I - 230.0  (sec) 

sp 

F - 444.822  (Nt) 

Cmax 

6 - 0.0  (rad) 

THRUST  VALUES  - 16850.0,  16850.0,  18471.5,  18471.5,  1621.5,  1621.5, 
0.0  (Nt) 

THRUST  TIMES  - 0.0,  0.200,  0.201,  1.650,  1.651,  13.350,  13.351  (sec) 
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No  Anomalies 

In  Simulation  1 the  turret  is  non-rotating,  there  is  no  tipoff,  no 
wind  and  the  rocket  is  "perfect,"  i.e.,  no  mass  imbalance  nor  thrust  mis- 
alignment. The  results  for  this  simulation  are  shown  in  Figs.  5-1  through 
5-17.  The  first  four  figures,  Figs.  5-1  through  5-8,  show  At^,  QNR  , QNR, 

Q,  Acz»  , Rjj^and  R as  functions  of  time.  These  figures  illustrate  the 

guidance  and  control  process.  As  the  rocket  leaves  the  launch  tube,  the 
force  due  to  gravity  causes  it  to  drop  and  hence  a pitch  up  command  is 
needed  (see  Fig.  5-2).  However,  when  the  rocket  leaves  the  launch  tube, 
guidance  is  not  possible  because  the  rocket's  aft  end  is  not  within  the 
FOV  of  the  IR  sensor.  From  the  time  of  sustainer  motor  ignition  to  the 
time  the  rocket's  aft  end  is  within  the  FOV  of  the  IR  sensor,  Q,Tr.  and  R„„ 

NRj,  NRj, 

or  zero  and  the  rocket  is  only  rate  stabilized.  Since  its  transverse 
angular  rates  at  the  beginning  of  the  flight  phase  were  zero  for  this  perfect 
launch  situation,  no  action  was  taken  by  the  control  system. 

For  this  simulation,  guidance  commands  were  possible  after  about  0.40 

seconds.  Due  to  the  rocket's  orientation  at  that  time,  lateral  control 

was  first  exerted,  then  longitudinal  control.  Note  that  At^  and  Atz  had 

maximum  magnitude  at  the  start  of  guided  flight,  but  not  for  the  total 

simulation.  Note  that  QKTD  (Fig.  5-2)  is  composed  of  "pulses"  for  which 

NK 

C 

the  maximum  values  vary  smoothly  with  time.  Note  also  that  the  amplitude 
of  the  oscillatory  curve  obtained  by  connecting  the  midpoints  of  the  ex- 
tremities of  adjoining  steps  has  the  appearance  of  a damped  oscillation. 


The  time  history  of  the  angular  rate  about  the  y^  axis,  Q^,  is  shown 
in  Fig.  5-3.  Note  the  rapid  response  in  Q,_.  Also,  of  interest  are  the 
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Figure  5-3.  Non-rolling  Angular  Rat 
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small  "bumps"  and  the  Q . curve.  These  are  due  to  the  application  of 

control  forces  during  lateral  control  phases,  since  a finite  control  region 

centered  on  <|>  + 6 = [ (2n+l)/2]7r,  for  example,  with  At^  positive,  results  in 

first  a small  negative  torque  about  that  axis  and  then  a positive  torque. 

The  net  effect  of  each  pair  of  positive  and  negative  torques  is  essentially 

zero.  We  note  that  if  there  were  no  aerodynamic  damping  of  transverse 

angular  rates,  the  tops  of  the  "pulses"  in  would  (except  for  the  bumps) 

be  flat  between  control  force  applications. 

The  angular  rate  about  the  y-axis  is  shown  in  Figure  5-4.  Because  6-0 

for  these  simulations,  the  control  axis  is  the  y-axis.  Therefore,  rapid 

changes  in  Q with  time  occur  when  control  of  either  a longitudinal  or  lateral 

motion  of  the  rocket  takes  place. 

The  results  for  Atz,  and  R are  similar  to  those  for  At  , 

Z c y c 

and  Qnr,  respectively.  However,  the  rates  R^R  and  R^R  are  generally  smaller 
!i ' c 

than  their  "longitudinal"  counterparts.  The  angular  rate  about  the  z-axis, 

R,  is  a much  smoother  function  of  time  than  Q,  since  it  results  from  the 

integrated  effect  of  control  torques  applied  about  the  y-axis. 

Figures  5-9  and  5-10  are  plots  of  the  angles  and  6x  as  functions  of 

time.  Since  the  wind  was  zero  for  this  simulation,  a = a and  8 = 0 and 

* x x 

hence,  these  are  the  time  histories  of  the  angle  of  attack  and  side  slip 
angle,  respectively,  of  the  rocket.  Note  that  a and  8 are  90°  out  of  phase, 
as  should  be  the  case. 


The  pitch  and  yaw  angles  of  the  rocket  as  functions  of  time  are  shown 


in  Figs.  5-11  and  5-12.  In  this  simulation  the  "target"  was  stationary  and 
could  be  considered  to  be  located  on  a line  of  sight  with  an  elevation  of 
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18.5°.  The  rocket's  pitch  angle  was  therefore  18.5°,  initially.  Due  to  the 
"gravity  drop"  as  the  rocket  left  the  launch  tube  with  zero  aerodynamic 
incidence  angle,  the  guidance  and  control  systems  operated  to  produce  an 
increase  in  0 for  a period  of  time  during  which  the  rocket  reacquired,  and 

i 

overshot  the  line  of  sight.  Long  period  oscillations  in  0,  with  smaller 
short  period  oscillations  superimposed  then  occurred.  The  yaw  angle  ip  was 
always  small  with  mean  value  approximately  zero.  Figure  5-13  is  included 
because  it  shows  clearly  the  maximum  and  minimum  values  of  0 and  which 
occurred. 

The  x -coordinate  of  the  center  of  mass  of  the  rocket  as  a function 

l J 

of  time  is  presented  in  Fig.  5-14.  It  is  essentially  a cubic  function  of 
time  during  the  portion  of  flight  simulated. 

Three  "views"  of  the  trajectory  of  the  center  of  mass  of  the  rocket 
as  seen  in  the  Gx„y  z system  are  given  in  Fig.  5-15.  Also  shown  in  Fig.  5-15 

w u u 

is  the  boundary  of  the  1°  FOV  of  the  IR  sensor.  This  is  the  smaller  of  two 
FOV's  used,  the  larger  being  4.5°.  We  note  that  although  the  transverse 
angles  are  not  very  smooth  functions  of  time,  the  trajectory  of  the  rocket's 
center  of  mass  is  very  smooth.  The  "head-on"  view  given  in  Fig.  5-15  shows 
that  the  center  of  mass  of  the  rocket  stays  within  2 meters  of  the  launcher 
to  target  LOS  during  the  course  of  the  simulation. 

Tipoff 

For  the  simulation  in  which  tipoff  from  the  launcher  was  included  and 
for  Simulations  3-10,  only  the  three  views  of  the  trajectory  are  given.  The 
three  views  of  the  trajectory  for  the  tipoff  simulation  are  shown  in  Fig. 5-16. 
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Lgure  5-15.  Trajectory  for  Simulation 


Lgure  5-16.  Trajectory  for  Simulation 
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The  main  effects  of  tipoff  are,  of  course,  that  the  rocket  has  in  initial 
transverse  angular  rate  and  its  center  of  mass  is  also  moving  transversely 
when  it  clears  the  launch  tube.  Upon  comparing  the  trajectory  for  the 
nominal  simulation  (Simulation  1)  with  that  given  in  Fig.  5-16,  we  note  little 
difference  in  the  two  sets  of  figures.  Hence,  by  itself,  tipoff  (at  least 
for  a distance  less  than  or  equal  to  0.15  mjis  not  significant. 

Mass  Imbalance 

Since  no  rocket  is  perfect,  each  may  be  expected  to  be  dynamically 
unbalanced.  Simulation  3 was  made  to  determine  whether  such  mass  imbalance 
would  have  a significant  effect  on  the  rocket’s  trajectory.  In  the 
simulation  code  the  angles  y ^ and  hj  (which  are  small  angles  of  rotation 
about  the  y-  and  z-axes,  respectively)  were  used  to  generate  time  varying 
products  of  inertia, 


I = -y,(I  -I  ) 
xy  3 yy  xx 


(5-la) 


I = y_(I  -I  ). 
xz  2 zz  xx 


(5-lb) 


For  yp=yy=  .0007,  the  three  views  of  the  trajectory  of  the  rocket's  center  of 
mass  shown  in  Fig.  5-17  were  produced  via  simulation.  This  particular 
source  of  error  above  does  not  appear  to  present  any  problem  as  far  as 
guidance  and  control  of  the  rocket  are  concerned. 

Thrust  Misalignment 

Figure  5-18  shows  three  views  of  the  trajectory  generated  by  including 
only  thrust  misalignment  as  an  error  source.  For  the  values,  ay  “ az  = .00035, 
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we  note  that  thrust  misalignment  is  also  not  a severe  problem.  That  is, 
although  there  are  recognizable  differences  in  the  trajectories  for 
Simulations  1 and  4,  these  are  not  large  enough  to  cause  concern. 

Cross  Wind 

Several  simulations  were  conducted  to  determine  the  relative  effects 


of  cross  winds  on  the  motion  of  the  rocket.  Simulation  5 is  representative 
of  the  results  obtained.  In  5,  a steady  cross  wind  of  6 meters/second  in 
the -Y -direct ion  was  included  in  the  input  for  an  otherwise  nominal  run.  By 
comparing  the  trajectory  shown  in  Fig.  5-19,  with  the  nominal  trajectory 
(Fig.  5-15)  small  differences  may  be  noted.  In  particular,  the  rocket 
(since  it  is  statically  stable)  turns  into  the  wind.  This  is  seen  clearly 
by  comparing  Figs.  5-15  and  5-19. 

Turret  Angular  Rates 

To  track  a moving  target,  the  turret  of  the  launch  must  rotate  about 
the  vertical  axis  z^.  The  rapidity  of  this  rotation  depends  on  the  speed 
of  the  target,  its  distance  from  the  launcher  and  its  direction  of  flight. 
The  Gx^ygZg  system  also  rotates  about  its  yG~axis  during  the  tracking 
process.  In  this  study,  we  have  neglected  this  rate  of  rotation  in  the 
launch  phase  simulations,  but  not  in  the  flight  phase  simulations, 
jr  For  Simulations  6 and  7,  the  target’s  motions  were  such  that  the  initial 

(for  flight  phase)  angular  rates  of  the  launcher  turret  about  the  zT-axis 

f 

were  +5°/sec  and  +10°/sec,  respectively.  The  target's  initial  position  for 

( 

the  two  cases  was  (1800,  0,  -600)  km  and  the  velocities  were  [302  157  0]T 

T 

m/sec  and  [130  314  0]  m/sec,  respectively.  The  rocket  was  fired,  in  both 
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cases,  from  the  right-hand  tube  of  the  launcher.  Thus,  in  these  simulations 
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the  FOV  of  the  IR  sensor  was  rotated  so  that  it  tended  to  "engulf"  the 
rocket.  This  resulted  in  the  entry  of  the  rocket  into  the  4.5°  FOV  more 
quickly  than  in  the  nominal  case  (Simulation  1).  Hence,  guidance  was 
begun  sooner  after  launch  in  Simulations  6 and  7 than  in  the  nominal  case. 

The  trajectories  for  Simulations  6 and  7 are  shown  in  Fig.  5-20  and  Fig.  5-21, 
respectively.  These  results  indicate  that  the  system  simulated  would  per- 
form adequately  even  for  the  10°/sec  turret  angular  rate  if  the  rocket  were 
fired  from  the  right-hand,  or  "leading"  tube. 

Results  for  negative  turret  angular  rates  about  the  zT~axis  are  shown 
in  Fig.  5-22  and  Fig.  5-23,  for  cases  in  which  the  initial  rates  were  -5°/sec 
and  -10°/sec,  respectively.  These  were  obtained  by  using  the  data  used  in 
Simulations  6 'hnd  7,  except  the  target  velocities.  These  velocities  were 
replaced  by  [301  -157  0]T  and  [130  -314  0]^  m/sec,  respectively,  for  Simu- 
lations 8 and  9.  The  rocket  was  acquired  by  the  IR  sensor  when  the  initial 
rate  was  approximately  -5°/sec,  since  its  aft  end  did  come  within  the  wider, 
4.5°  FOV  of  the  sensor  (see  Fig.  5-22).  Also,  the  guidance  and  control 
system  quickly  brought  the  rocket  within  the  narrower,  1°  FOV  of  the  sensor 
and  it  remained  therein.  However,  in  the  -10®/sec  case,  the  IR  sensor's 
FOV  was  rotated  away  from  the  rocket's  initial  direction  of  flight  so 
rapidly  that,  although  the  rocket  entered  the  4.5°  FOV  momentarily,  it  did 
not  remain  therein.  Hence,  guidance  and  control  of  the  rocket  was  lost. 
Combination  of  Anomalies 

Simulation  10  was  made  with  a combination  of  anomalies  as  indicated 
in  Table  5-1.  As  can  be  seen  in  Fig.  5-24,  the  rocket  entered  the  4.5°  FOV 
and  remained  within  the  1°  FOV  after  reaching  a range  of  approximately  310  m. 


0.0 


gure  5-23.  Trajectory  for  Simulation 
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Figure  5-24-  Trajectory  for  Simulation  10. 


6.  SUMMARY,  CONCLUSIONS  AND  SUGGESTIONS 
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6.1  Summary 

In  this  report  we  have  described  the  physical  and  mathematical  models 
which  form  the  basis  for  digital  simulation  codes  which  may  be  used  in 
studying  the  motion  of  a short-range  air  defense  system  rocket  during 
launch  and  prior  to  its  acquisition  by  a launcher-based  radar.  Listings 
of  the  FORTRAN  IV  digital  computer  codes  are  also  included.  The  codes  were 
used  to  simulate  a typical  system  and  results  of  these  simulations  have 
been  presented  and  discussed. 

6.2  Conclusions 

The  simulation  codes  are  flexible  enough  that  different  systems  may  be 
simulated  by  changing  data  and/or  subroutines.  They  are  not  as  complex  as 
some  codes  of  a similar  nature,  since  simulation  of  the  operation  of  all 
the  various  subsystems  is  not  attempted.  However,  many  anomalies  are 
modeled  and  the  computer  time  required  for  each  run  is  reasonable  (around 
2 .5  minutes  on  an  IBM  370/158  in  conjunction  with  the  FORTRAN  H compiler. 

In  regard  to  the  particular  system  simulated,  it  was  found  that  the 
presence  of  anomalies  such  as  cross-winds,  mass  imbalance,  thrust  misalign- 
ment, etc.,  of  expected  magnitudes  did  not  prevent  adequate  operation  of 
the  system.  The  only  factor  which  seems  to  be  of  major  importance  is  the 
rate  of  rotation  of  the  turret,  and  this  is  of  vital  importance  only  when 
the  rate  is  large  in  magnitude  (say,  10°/sec)  and  the  rocket  is  fired  from 
a "lagging”  tube.  This  would  appear  to  preclude  the  use  of  both  rockets 
against  two  fast  moving  targets  which  appear  successively.  That  is,  if  a 
target  moving  such  that  a large  turret  rate  is  necessary  appears,  if  a rocket 
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is  fired  from  the  leading  tube  and  the  rocket  guided  to  impact  and,  if  before 
the  time  required  for  reloading  has  elapsed,  another  such  target  appears, 
based  on  the  simulations  made  in  this  study,  the  second  target  cannot  be 
successfully  attacked. 

i 

6.3  Suggestions 

On  the  basis  of  experience  gained  during  the  course  of  this  study,  it  . 

is  suggested  that  additional  attention  be  given  to  determining  the  effects 

of  the  following  factors:  | 

1.  Induced  motions  of  the  rocket  due  to  launcher  motions  j 

other  than  the  angular  rate  of  the  turret  about  a J 

vertical  axis.  1 

2.  Unfolding  of  the  wings  of  the  rocket.  j 

It  is  also  suggested  that  the  simulation  codes  presented  herein  be  j 

modified  and/or  extended  to  include  mathematical  models  of  the  above 

factors.  j 
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APPENDIX  A 


LISTING  OF  MISSILE  SIMULATION  (MISSIM)  PROGRAM 
On  the  following  pages  is  a listing  of  the  missile  simulation  (MISSIM) 
program  used  for  the  flight  phase  analysis  in  the  report.  The  computer 
code  developed  is  compatible  with  the  FORTRAN  G,  FORTRAN  H,  or  WATFIV 
compilers  used  in  conjunction  with  the  IBM  370/158  digital  computer  avail- 
able at  Auburn  University.  All  computations  in  the  program  are  done  in 
the  double  precision  mode. 

All  input  into  the  program  is  formatted  with  a generalized  (G)  field 
descriptor  and  is  accomplished  with  punched  cards.  Each  card  contains 
ten  data  fields  which  are  eight  columns  wide.  The  generalized  descriptor 
was  chosen  because  it  allows  one  to  input  logical,  integer  or  real  (external 
fixed  point  or  floating  point)  variables  using  the  same  format.  The  entire 
input  necessary  to  use  MISSIM  is  read  in  the  main  program,  and  is  explained 


below. 


NVAL  The  number  of  thrust-time  values  tabulated  for  the  missile 

(right-hand  justified). 

TVECTR  Vector  (dimensioned  NVAL)  of  thrust  magnitude  times  (sec) . 

FVECTR  Vector  (dimensioned  NVAL)  of  thrust  magnitude  values  (Nt). 

XL  Length  of  the  missile  (m) . 

S Reference  area  for  aerodynamics  (m2). 

D Reference  length  for  aerodynamics  (m), 

RP  Outer  diameter  of  the  booster  motor  propellant  grain  (m). 


Initial  inner  diameter  of  the  booster  motor  propellant 
grain  (m). 
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RHOP 

PGL 

ISP 

GRAV 

RHO 

TEMP 

GAMMA 

RAIR 

XP1.XP2 

SLAT 

RT 

TIMEV 

VXCG 

VIX 

VIY 

VIZ 

IXY,  IXZ, 
IYZ 

FCMAX 

DELTA 

RTUR 


Booster  motor  propellant  density  (kg/  3) . 

Booster  motor  propellant  grain  length  (m) . 

Booster  motor  special  impulse  (sec) . 

Gravity  magnitude  (m/sec2) . 

Atmospheric  ambient  density  (kg/m3). 

Atmospheric  ambient  temperature  (°K) . 

Ratio  of  specific  heats  for  air. 

Gas  Constant  for  air  (m2/s2-°X) . 

Coordinates  of  ends  of  booster  motor  propellant  grain  (m) . 

Distance  from  point  A to  aft  end  of  missile  (m) . 

Vector  from  point  A to  a point  on  the  line  of  action  of 
the  sustainer  motor  thrust  (m). 

Five-element  vector  of  times  for  CG  location  and  moments  of 
inertia  values  (sec) . 

Five-element  vector  of  CG  locations  (m). 

Five-element  vector  of  spin  moment  of  inertia  values  (kg-m2). 

Five-element  vector  of  pitch  moment  of  inertia  values  (kg-m2). 

Five-element  vector  of  yaw  moment  of  inertia  values  (kg-m2) . 

Values  of  products  of  inertia.  (Used  if  mass  imbalance 
quantities  are  not  specif ied)  (kg-m2 ) . 

Magnitude  of  control  force  (Nt) . 

Angle  that  the  control  plane  makes  with  the  z-axis  (see 
Fig.  3-4) (rad). 

Distance  from  axis  of  rotation  of  launcher  turret  to  origin 
of  goniometer  system  (m)  . 
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HANT 


Height  of  radar  antenna  from  origin  of  inertial  system  On). 
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K 


TBEAM  Half-angle  of  infrared  seeker's  field  of  view  (rad). 
TMAX  Maximum  allowed  thrust  deflector  dwell  time  (sec) . 


TDELAY 

MAXRAT 

TBOOST 

TBC 

EULER 

PRTSUP 

PLONLY 

H 

NEQ 

NUMB 

NUMPPT 

UG,VG,WG 

IF 

TO 

ALY.ALZ 

MU2,MU3 


Delay  time  from  command  of  control  force  until  control 
force  magnitude  equals  FCMAX  (sec). 

Macimum  allowed  commanded  angular  rates  (rad /sec). 

Time  of  booster  motor  burnout  (sec) . 

Time  of  sustainer  motor  ignition  (equals  earliest  time  that 
control  is  possible) (sec) . 

Logical  variable  defining  type  of  integration  routine 
desired.  If  Euler  is  true,  integration  is  Euler;  otherwise, 
fourth-order  Runge-Kutta  integration  is  used. 

Logical  variable  indicating  whether  excess  output  is  to  be 
suppressed. 

Logical  variable  indicating  whether  output  is  to  be  printed 
only  on  the  dataset. 

Integration  step  size  (sec) . 

Number  of  dependent  variables  to  be  integrated  (right- 
hand  justified). 

Number  of  integration  steps  before  output  is  printed  (right- 
hand  justified). 

Number  of  integration  steps  before  output  is  stored  for 
Dataset  storage  (right-hand  justified). 

Components  of  wind  velocity  in  negative  direction  of  inertial 
axes  (m/sec). 

Time  at  end  of  simulation  (sec)  . 

Time  at  beginning  of  simulation  (sec). 

Euler  angles  which  orient  the  thrust  vector  (rad ) . 

Euler  angles  which  define  principal  axis  rotations.  Used  to 
compute  dynamic  mass  imbalance  quantities  (rad). 


X 


Vector  of  initial  values  of  the  dependent  variables 
(dimensioned  NEQ) . 
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XTV 

Initial  position  vector  of  target  in 
system  (ra) . 

inertial  reference 

VTV 

Velocity  vector  of  target  in  inertial 
(m/sec) . 

reference  system 

GAIN1 

Control  gain  K1  (sec2). 

GAIN2 

Control  gain  K2  (sec/m). 

GAIN3 

Control  gain  K3  (sec2/m)« 

TITLE 

Eighty-character  alphanumeric  vector 
title  of  the  simulation. 

for  the  identifying 

Several  additional  vectors  of  input  data  used  in  subroutine  AERO  to  de- 
termine the  aerodynamic  force  and  moment  on  the  missile  are  described  in 
Appendix  C. 

Output  from  the  program  can  be  obtained  in  three  forms.  Which  form 
is  obtained  depends  on  input  values  of  PRTSUP  and  PLONLY  (explained  above). 
In  the  normal  form  (PRTSUP  = FALSE),  much  output  is  printed  for  checking 
purposes.  In  the  suppressed  form  (PRTSUP  = TRUE  and  PLONLY  = FALSE),  only 
the  time  and  the  state  variables  are  printed.  In  the  third  form  (PRTSUP  = 
TRUE  and  PLONLY  = TRUE)  output  is  stored  on  DATASET  only,  and  no  printed 
output  is  obtained. 

Comment  cards  have  been  provided  liberally  throughout  the  program  to 
indicate  the  most  important  calculations  and  to  provide  the  user  with 


additional  information. 
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NOW  THAT  THE  K'S  HAVE  BEEN  CALCULATED,  COMPUTE  THE 
VALUE  Of  DELI  AX  TROtt  THE  WE10H1ED  AVERAGE  FORMULA. 
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DI.MEN.SICN  A ( »,  J)  , ;t  ( ».  1)  .SPill.CPHI) 

DJUdLE  PkECIBION  A,  D,UET, D, DADB, AuJD  IMPLICIT  REAL*D  (A-H.O-2) 

DO  5 I = 1,J  DIMENSION  C ( J , J)  , CT  (3,  3)  , V ( 2k)  , TM1  ( 3 , i)  , TriE  (3.  J) 
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LISTING  OF  SIMPLIFIED  MISSILE  LAUNCHER 
(MISSLNCH)  SIMULATION  PROGRAM. 

On  the  following  pages  is  a listing  of  the  missile  launcher  (MISSLNCH) 
computer  program.  MISSLNCH  is  used  to  simulate  the  missile  dynamics  during 
the  launch  phase  of  the  analysis  in  order  to  obtain  initial  conditions  for 
the  flight  phase  simulation  (MISSIM)  program.  The  computer  code  developed 
is  compatible  with  the  FORTRAN  G,  FORTRAN  H,  and  WATFIV  compilers  used  in 
conjunction  with  the  IBM  370/158  digital  computer  available  at  Auburn  Univer- 
sity. All  computations  in  the  program  use  double  precision  arithmetic. 

All  input  into  the  program  is  accomplished  by  punched  cards.  Each 
card  contains  ten  data  fields  which  are  eight  columns  wide.  A generalized 
(G)  field  descriptor  is  used  in  the  format  for  all  read  statements.  This 
format  was  chosen  for  simplicity,  since  it  allows  integer,  logical,  or  real 
(external  fixed  point  or  floating  point)  variables  to  be  read  with  the  same 
format.  The  entire  input  necessary  to  use  MISSLNCH  is  read  in  the  main 
program  and  is  explained  below. 


IS 

Spin-axis  moment  of  inertia  (kg-m2). 

IT 

Transverse-axis  moment  of  inertia  (kg-m2). 

THETAT 

Elevation  (or  pitch)  angle  of  launch  tube  (rad). 

OMEGAT 

Steady  angular  velocity  of  launch  turret  measured 
vertical  axis  (rad/sec). 

about 

XA , YA 

X and  Y components,  respectively  (measured  in  the  turret 
reference  system),  of  any  point  fixed  within  the  launch 
tube  relative  to  a point  on  the  turret  axis  of  rotation 
(m) . 
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DIA  Missile  largest  diameter  (m). 

MASSR  Missile  mass  (kg). 

GRAV  Magnitude  of  gravity  (m/sec2). 

FT  Steady  thrust  magnitude  which  acts  on  missile  (Nt). 

TX  Spin  torque  applied  to  missile  during  launch  (Nt-m) . 

XETO  The  distance  from  any  point  fixed  within  the  launch  tube 

to  the  missile  center  of  mass  when  tipoff  has  ended  (m) . 

TODIST  The  distance  that  the  missile  travels  during  the  tipoff 
interval  (m) . 


XGT.ZGT 
X 

PRMT(l) 

PRMT(2) 

PRMT  (3) 

PRMT(4) 

ERROR 

Since  this  computer  program  is  to  be  used  only  during  the  short  time 
during  the  launch  phase,  the  missile  physical  model  is  considerably  more 
simple  than  that  used  in  the  flight  phase.  During  the  launch  phase,  the 
missile  is  modeled  as  a rigid,  constant  mass  body  with  constant  moments  of 
inertia,  which  is  acted  upon  by  a constant  thrust  and  initial  spin  torque. 

In  addition,  the  MISSLNCH  program  which  follows  is  also  a simple  mathemathcal 
model  of  the  missile/launcher  system,  and  does  not  contain  the  level  of 


The  X and  Z components,  respectively  (measured  in  the  tube 
reference  system),  of  point  G,  the  center  of  the  infrared 
angle  sensor  (m) . 

The  vector  of  initial  conditions  of  the  state  variables. 

Time  at  the  beginning  of  the  simulation  (sec) . 

Time  at  the  end  of  the  simulation  (sec). 

The  maximum  allowed  integration  step  size  (sec) . 

An  upper  error  round  to  test  for  convergence  in  the  inte- 
gration algorithm. 

A vector  containing  error  weights  of  the  state  variables. 

Is  used  to  test  for  convergence  in  the  integration  algorithm 
(The  sum  of  the  components  must  equal  1) . 
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generality,  nor  the  inclusion  of  large  numbers  of  comment  cards  to  explain 
the  sections  of  the  program  that  was  evident  in  the  MISSIM  program  of 
Appendix  A.  If  a more  complex  missile/launcher  model  is  desired,  the 
reader  is  directed  to  Ref.  2 for  a listing  of  *Tie  more  complicated  math- 
ematical model  used  during  that  study. 

Some  subroutines  used  in  MISSLNCH  are  also  used  in  MISSIM.  In  order 
to  avoid  repetition,  their  listings  are  presented  only  in  Appendix  A.  The 
subroutines  are  TILDE,  TRANSP,  MATINV,  MATMPY , MATXV , CXMAT,  CXVEC,  VAC, 
VSUB,  MAD  and  MSUB. 
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APPENDIX  C 


AERODYNAMIC  COEFFICIENTS 


In  this  appendix,  aerodynamic  data  for  a typical  SHORADS  rocket,  the 

motion  of  which  was  simulated  in  this  study,  is  presented  and  the  way  in 

which  this- data  is  used  in  the  flight  simulation  code  is  described.  This 

data  was  provided  by  Mr.  John  Howerton,  U.S.  Army  Research  and  Development 

Command,  Redstone  Arsenal,  Alabama. 

C.,  and  C 
N y 

The  coefficients  CN  and  C^  are  obtained  from  C^  for  which  values  are 

a 

given  in  Table  C-la  for  twenty-seven  Mach  number  values  and  two  incidence  angle 
(i^  = /i2+p2)  ranges,  i^  4°  and  iT  >_  4°.  During  the  simulation,  when  values 
of  CN  and  Cy  are  needed,  the  incidence  angle  is  computed  and  the  proper  one- 
dimensional array  of  stored  C^  values  is  chosen.  Corresponding  Mach 

a 

numbers  which  bracket  the  flight  Mach  number  are  determined  next  and  the 
value  of  C^  used  in  computing  C^  and  C^  is  determined  by  linear  interpola- 
tion. For  the  rocket  considered  herein,  C„  and  C for  i„  > 4°  are  computed 

N y T 

using  C^  = C^  + CN  a and  C = - CN  - C^  6,  where  CN  is  a function  of 
o a * o a o 

Mach  number  as  shown  by  Table  C-la. 


Values  of  the  axial  force  coefficient ,C,  , as  a function  of  Mach  number 

A 

o 

are  given  in  Table  C-la  also.  A linear  interpolation  subroutine  contained 

in  the  simulation  code  is  used  to  determine  the  proper  value  of  C^  . 

o 


TABLE  C-l.  AERODYNAMIC  CHARACTERISTICS. 


M 

CN 

a 

it- ° 

CN 

a 

1t>4° 

CN 

o 

it>4° 

CA 

0 

0.0 

15.0 

26.0 

-0.8 

0.44 

0.1 

15.0 

26.0 

-0.8 

0.44 

0.2 

15.0 

26.0 

-0.8 

0.433 

0.3 

15.6 

26.1 

-0.79 

0.430 

0.4 

15.0 

26.1 

-0.78 

0.42 

0.5 

15.5 

26.4 

-0.74 

0.416 

0.6 

16.0 

26.5 

-0.70 

0.41 

0.7 

16.25 

26.7 

-0.71 

0.41 

0.8 

16.2 

26.5 

-0.72 

0.41 

0.85 

16.25 

26.7 

-0.725 

0.415 

0.90 

16.7 

27.0 

-0.74 

0.435 

0.925 

17.1 

27.5 

-0.74 

0.46 

0.95 

17.5 

27.6 

-0.74 

0.50 

0.975 

17.7 

27.8 

-0.725 

0.55 

1.0 

18.0 

28.0 

-0.715 

0.605 

1.025 

18.2 

28.3 

-0.71 

0.645 

1.05 

18.5 

28.5 

-0.705 

0.665 

1.075 

18.6 

28.6 

-0.70 

0.685 

1.1 

19.0 

28.7 

-0.70 

0.69 

1.15 

19.1 

28.9 

-0.70 

0.68 

1.2 

18.9 

28.5 

-0.70 

0.66 

1.25 

18.8 

28.0 

-0.70 

0.605 

1.3 

18.6 

27.5 

-0.695 

0.57 

1.4 

18.5 

27.0 

-0.695 

0.55 

1.5 

18.0 

26.2 

-0.625 

0.54 

1.6 

18.0 

25.5 

-6.550 

0.53 

1.7 

18.0 

25.0 

-0.480 

0.53 

(a) 
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TABI.K  C~  1 (CONT) 


10P 


The  roll  damping  coefficient,  , is  given  in  Table  C-lb  as  a function 

P 

of  incidence  angle,  i^,  for  Mach  numbers  less  than  and  greater  than  unity. 

In  the  simulation  code  linear  interpolation  is  used  to  find  the  proper 

P 

value . 


The  pitch  and  yaw  rate  damping  coefficients, C and  C , respectively, 

mq  nr 

are  equal,  since  the  rocket  is  geometrically  symmetric  with  respect  to  the 

x-axis.  For  canard  (destabilizing  fins)  undeployed  and  deployed,  values 

of  the  coefficient  C are  given  in  Table  C-lc  along  with  corresponding 

values  of  l-x  , the  distance  of  the  center  of  mass  of  the  rocket  from  its 
eg’ 

nose.  The  distance  II  is  computed  in  the  simulation  code  during  each 
evaluation  of  the  derivatives  of  the  state  variables  and  the  linear  inter- 
polation subroutine  is  used  (as  in  other  calculation  discussed  above)  to 

find  the  proper  value  of  C = C 

m n 

q r 

C and  C 
m n. 

a 8 

For  the  rocket  simulated  herein,  C = - C and  C is  a function  of 

n.  m m 

p a a 

the  center  of  pressure  location  x^  (measured  from  the  aft  end  of  the 
rocket)  which  in  turn  is  a function  of  Mach  number,  aerodynamic  incidence 


angle  and  B^,  the  percent  of  the  total  area  of  the  canards  (destabilizing 

fins)  which  is  exposed  to  the  flow  field  surrounding  the  rocket.  In  the 

simulation  code,  C ■ C„  (x*  - x ) and  x*  ■ f. (B  ,M)  + f»(i_).  Here 
* m N cp  eg  cp  1 r’  2 T ’ 

a a 

Br  is  dependent  upon  £ , the  distance  the  deployment  piston  has  traveled. 

The  distance  l is  a linear  function  of  time  which  is  zero  at  0.8  seconds 
r 
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after  booster  motor  Ignition  and  has  its  maximum  value  at  1.7  seconds  after 


booster  motor  ignition. 


To  find  C , the  value  of  l is  first  determined.  Then  B is  determined 
m r r 

a 

by  linearly  interpolating  using  the  values  in  Table  C-ld.  Next,  f^(Br>M) 
is  found  by  using  a double  linear  interpolation  scheme  and  the  data  of 


' • Table  C-le.  The  value  of  i^,  and  the  corresponding  value  for  are  then 

found  from  the  data  listed  in  Table  C-lf.  These  values  of  f^  and  {2  are 

► used  to  find  x*  and  the  current  value  of  x is  used  along  with  the  x' 

cp  eg  cp 

value  and  the  appropriate  value  of  CN  to  compute  CM  . 

a a 


